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ABSTRACT 

We predict near-infrared luminosity functions of young (5 Myr to 1 Gyr) star 
clusters by combining evolutionary models of very low-mass (1 Mj to 0.15 Mq) 
dwarfs with empirical bolometric corrections. We identify several characteristic 
features in our results. These can be attributed to three causes: (1) deuterium 
burning in the most massive substellar objects; (2) methane absorption in bodies 
with Tg// less than 1300 K, the temperature of the L/T transition; and (3) the 
formation of dust clouds and the rainout of dust at roughly the same effective tem- 
perature as methane formation. Accurate reconstruction of the substellar mass 
function from luminosity function observations requires that these phenomena 
are taken into account. At present, few observational studies extend to sufficient 
sensitivities to allow detection of these effects. However, the luminosity function 
of the young open cluster IC 2391 shows a clear peak at Mj ~ 14 which we 
attribute to the result of deuterium burning in substellar objects. The location 
of this feature is a strong function of age, and we estimate an age of 35 Myr for 
IC 2391. This is significantly younger than the 53 Myr derived from the location 
of the lithium depletion boundary but agrees with the main sequence turnoff age. 
We consider the implications of this result and our multi-band luminosity func- 
tions for future observational studies. All predicted luminosity function features 
are, or will be, accessible to observations using new wide-field IR imagers and 
the Space Infrared Telescope Facility. 

Subject headings: stars: evolution — stars: low-mass, brown dwarfs — stars: 
luminosity function, mass function 
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1. Introduction 

The determination of the stellar Initial Mass Function (IMF) has been a prominent goal 
of astronomy for over 50 years. The IMF, ^(M), describes the number of stars born per 
unit mass per unit volume. The pioneering study by Salpeter (1955) was limited to stars 
with masses greater than 0.3 Mq. Salpeter found that the results were well represented 
by a power law, *(M) oc M"", with a = 2.35. Two decades later. Miller & Scalo (1979) 
used improved observations of lower luminosity stars to show that the IMF deviates from 
a simple power law below 1 Mq. They derived a better fit with a log- normal distribution, 
^(M) cx exp (i^5(i^^^i(^)2^ ^ith Mo ~ 0.15 Mq, a - 0.6, and a distribution close to 
the Salpeter power-law at masses exceeding 1 Mq. The change in slope at lower masses 
is evident both in field-star surveys and in observations of young open clusters (e.g., Scalo 
(1998)). More recent studies (Reid & Gizis 1997; Kroupa 2001; Reid, Gizis, & Hawley 2003) 
favor modehng the field-star IMF using multiple power-laws, with a ~ 2.3 at M > 1 Mq 
and a ~ 1 at lower masses. 

Extending coverage to masses close to and below the hydrogen-burning limit has only 
become possible within the last decade, as improvements in detector and computational 
technologies have permitted the first high-sensitivity, wide-field near-infrared sky surveys. 
Analysis of those surveys has led to the discovery of a host of very low-mass ultracool dwarfs 
(Kirkpatrick et al. 2000; Burgasser et al. 2002; Geballe et al. 2002) and the first estimate of 
the brown dwarf mass function in the field (Reid et al. 1999). 

There are a number of comphcations in deriving ^{M) for substellar mass objects. 
Mass is generally not observed directly, but is estimated from the measured luminosity. In 
the case of main-sequence stars this calculation is relatively straightforward, since hydrogen 
fusion leads to a well defined mass-luminosity relation. Brown dwarfs, however, lack a stable, 
long-term energy source, and, as a result, evolve rapidly to lower temperatures and lower 
luminosities. Thus, the substellar luminosity function depends strongly on both mass and 
age. Moreover, the presence of broad spectral absorption features stemming from the onset 
of molecule and dust cloud formation can lead to significant changes in broadband absolute 
magnitudes over small ranges of effective temperature. Finally, short periods of deuterium 
fusion occur in objects with masses greater than 13 Mj, leading to further transient effects 
in the luminosity function. 

We invert the problem. Starting with an assumed mass function and age distribution, 
we combine theoretical predictions of the evolution of very low-mass dwarfs from the models 
computed by Burrows et al. (2001) with empirical estimates of bolometric corrections from 
Dahn et al. (2002) and Reid & Cruz (2002) to predict luminosity functions for young (5 Myr 
to 1 Gyr) clusters. Young clusters are convenient for mass function studies because they 
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offer a large range of masses with similar ages and distances. Muench, Lada, & Lada (2000) 
used the Baraffc ct al. (1998) models to perform a similar analysis, to masses as low as 0.02 
Mq, and probed a complementary age range (< 10 Myr) to the current work. 

The outline of this paper is as follows: §2 provides a summary of the input evolutionary 
models and bolometric corrections for our analysis; §3 describes the main features of the 
predicted luminosity functions. In §4, we compare those results against observations of 
several young clusters, and consider the potential offered by future observations. 

2. Modeling the Substellar Luminosity Function 

2.1. SubstellEir Evolutionary Models 

Most current studies of very low-mass dwarfs are based on one of two sets of theoret- 
ical models. Baraffe et al. (1998) combine interior models with detailed model atmosphere 
calculations, but cover only a limited range of mass and effective temperature (essentially 
corresponding to a lower spectral type limit of L); Chabrier, Baraffe, AUard, & Hauschildt 
(2000) update these models to cover a lower range of masses, 0.1 Mq to 0.01 Mq (spectral 
types M, L, and T). In contrast. Burrows et al. (2001) present a more extensive series of 
models, spanning masses between 1 Mj and 0.15 Mq (1 Mq = 1047 Mj) and ages in the 
range 10^ to 10^° years, that do not include model atmospheres, but predict bolometric 
luminosities and effective temperatures as a function of time. Since our aim is to include 
coverage of very low-mass stars through low-temperature brown dwarfs, we base our analysis 
on the Burrows et al. (2001) models. 

Figure 1 shows the time variation in effective temperature predicted by the Burrows 
models. These models provide a series of reference points for specific masses at particular 
ages. To achieve the numerical resolution needed for our study, we interpolate among the 
predicted values of bolometric luminosity and effective temperature using a bicubic spline 
interpolation as a function of mass and age. This interpolation steps over fine increments 
in both mass (Am = IxIO'^Mq) and age (Alog(Age(Gyr)) = 8.602x10-''). We do not 
attempt to model clusters with ages of less than 2-3 Myr, since the model calculations are 
unreliable at these ages (Burrows et al. 2001). 

Our goal is to compare the model predictions against observations of near- and mid- 
infrared (IJHKM) luminosity functions of young stellar clusters. As a first step, bolometric 
luminosity functions are generated. This process requires assumptions concerning both the 
IMF and the distribution in age. We parameterize the IMF as a two-segment power-law, 
{ai,a2). The higher- mass segment includes the range 0.08 Mq-0.15 Mq, with ai — 1.05 
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held fixed (following Reid & Gizis (1997)). At lower masses, from 1 Mj to 0.08 Mq, is 
allowed to vary as a free parameter. The two segments of the power law are connected at 
the joining mass, m\2- For the set of models presented here mi2 is fixed at 0.08 Mq (see 
§4.1). 

The age distribution assumption is simplified because a single age is assigned to a 
cluster. However, star formation is not an instantaneous event, but rather a continuous 
process that occurs over a finite time interval. We assume that the ages for stars/brown 
dwarfs in each cluster uniformly span a range corresponding to ±10% of the assigned cluster 
age. Although this assumption produces an unrcalistically large age dispersion for older 
clusters, the luminosity variations over the 10% age interval should not be large enough to 
affect our predictions. 

Combining the IMF and the assumed age range gives a continuous distribution of objects 
as a function of mass and age. Each model mass-age distribution is normalized such that the 
space density of 0.1 M© objects is equal to that derived in Reid & Gizis (1997). This arbitrary 
normalization is used throughout this work unless otherwise stated. Lastly, we use the 
Burrows et al. (2001) models to assign values for the appropriate bolometric luminosity and 
effective temperature at each point in the mass-age distribution. The luminosity function is 
derived from this distribution by summing the relative number of objects in each interval of 
luminosity. 

2.2. Bolometric Corrections 

In order to compare model luminosity functions to observed data, our derived bolo- 
metric luminosity functions must be converted into broad-band luminosity functions. This 
conversion requires the use of bolometric corrections (BCs). Empirically derived BCs are 
more reliable than colors generated by model atmospheres, e.g., Baraffe et al. (1998) (for L 
and M dwarfs) or Burrows et al. (2001). Model BCs arc only as good as the current knowl- 
edge of line opacities of the atoms, ions, and molecules present in brown dwarf atmospheres, 
whereas empirical BCs have no such model dependencies. 

There are two sets of empirical near-infrared very low-mass star and substellar BCs 
available (Leggett et al. 2002; Dahn et al. 2002). The temperature scales derived by these two 
groups differ slightly, with the Dahn effective temperatures closer to those for the Burrows 
models. We therefore use the J-band BCs from Dahn et al. (2002) and derive Ji^i^-band 
BCs from the Dahn et al. (2002) IJEK colors. However, the Dahn et al. (2002) BCs lack 
points in the T dwarf regime, T < 1300K, so we include the data on Gl 229B and Gl 570D 
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from Leggett et al. (2002) in our analysis. For the M-band we use BCs from Table 2 of 
Rcid & Cruz (2002), plus spectrophotometric data on Gl 229B from Noll, Geballe, & Marley 
(1997). 

A combination of polynomials and line segments was used to create a fit of BC as a 
function of effective temperature in each bandpass (Figure 2). T dwarf data (the coolest 
one M-band and two i^-band points) are joined by straight line segments to the rest of 
the fit to accommodate the sharpness of the L/T transition. The effect at K-hand of the 
rapid onset of methane formation is seen at the L/T transition temperature of 1300 K 
(Figure 2, top panel). This effect leads to a sharp reduction in the i^-band luminosity over 
a relatively small range in effective temperature, with a correspondingly significant effect on 
the K-hsoid luminosity function, as is discussed in more detail in §3.2. Below the last data 
points available we extrapolate a fiat fine segment. This reflects our ignorance of the BCs 
below ~900 K. However, it is possible that additional structure exists at faint magnitudes 
in the luminosity function. We use our derived IJHKM BC fits to derive broadband 
magnitudes for each distribution of objects in our simulations and construct the IJHKM- 
band luminosity functions in the same manner as the bolometric luminosity functions. 



3. Model Results: Cluster Luminosity Function Morphology 

3.1. The Bolometric Luminosity Function 

Through application of the methods described in §2, we construct a library of model 
luminosity functions based on objects in the Burrows et al. (2001) mass range with ages 
between 5 Myr and 1 Gyr. We show in Figure 3 an example of the general results: the 
evolution of the bolometric luminosity function of a cluster from 10 Myr to 1 Gyr for a 
model with power law index 02 = 1.0. We identify in these luminosity functions significant 
features that are discussed below. 

At 10 Myr, the cluster exhibits a strong peak (A) at bright magnitudes, a broad trough 
(B) centered on M^o/ ~ 14, and a broader peak (C) at fainter magnitudes. Peak A is the 
low-mass tail of the main sequence (M > 0.075 Mq) together with brown dwarfs sufficiently 
massive to burn deuterium (13 Mj < M < 0.075 Mq). Note that the sharp cutoff at bright 
magnitudes is due to the high mass limit of the Burrows models, 0.15 M0, and is therefore 
unphysical. The presence of deuterium fusion in higher-mass brown dwarfs maintains their 
surface temperatures and luminosities at relatively high levels (sec Figure 1). In contrast, the 
non-dcutcrium burning brown dwarfs evolve rapidly to lower temperatures and luminosities. 
This causes the two populations to spread apart, forming a trough (B). Peak C represents 
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the lowest mass objects (M < 13 Mj), which cool continuously after their formation. 

The 100 Myr cluster luminosity function has a feature which is not present in the 10 
Myr luminosity function: a sharp spike (D) between peaks A and C. By 100 Myr, nearly all 
brown dwarf deuterium burning has ceased, and these higher-mass brown dwarfs join their 
lower-mass cousins in rapid cooling. This leads to the formation of the second peak/trough 
feature (D and E) as brown dwarfs separate from hydrogen-burning stars. To verify that 
peak D objects correspond to higher-mass brown dwarfs, we segregate the 100 Myr sample 
by mass. Figure 4 confirms that peak D includes only brown dwarfs that have finished their 
deuterium burning phase. 

Finally, the bottom panel in Figure 3 shows our prediction for a 1 Gyr cluster. The 
evolution of peak A has halted, since nearly all stars are stabilized on the main sequence 
and all brown dwarfs that burned deuterium have separated from this group. Peak D has 
broadened since brown dwarfs cool at mass-dependent rates. A new trough (F) lies in the 
midst of this peak, reflecting the onset of dust cloud formation and dust rainout in the 
Burrows models (see §3.2). Note that peak C spans roughly the same magnitude range at 
each age (all panels of Figure 3). In general, the cooling rates of substellar objects are mass 
dependent, but the cooling rates of the lowest mass (< 13 Mj) brown dwarfs are similar to 
each other between 5 Myr and 1 Gyr. This leads to the uniform behavior of peak C members 
over time. Wc do not include dynamical evolution effects in our models. While these effects 
arc not hugely important for the younger clusters considered here the effects on old clusters, 
~1 Gyr, may be important. 

Overall, the bolometric cluster luminosity function exhibits complicated time dependent 
behavior that depends on both the interior and atmospheric evolution of brown dwarfs. 
However, the features identified in Figure 3 are dependent mostly on the evolution of the 
total emergent luminosity (i.e. interior evolution). Therefore, since the interior physics of 
these models is relatively well understood, we believe that the locations of the major features, 
peak D and troughs B and E, are defined robustly. We note that the normahzation of the 
model luminosity functions ultimately is arbitrary and that only the relative heights and 
locations of these features are important. 

3.2. Broad-band Luminosity Functions 

Broad-band luminosity functions are constructed using the BCs discussed in §2.2. Figure 
5 displays the evolution of /- and K-hsaid cluster luminosity functions from 10 Myr to 
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900 Myr^. The morphological features discussed in §3.1 and labeled in Figure 3 are also 
marked in Figure 5. The gross morphological features of both bands are similar to those of 
the bolometric luminosity function. Peak A lies at the brightest magnitudes modeled here; 
the sharp cutoff is due to the upper mass limit of the Burrows models. Peak D (representing 
brown dwarfs that burned deuterium) forms and separates from the (hydrogen-burning) 
main sequence at ages of 10-20 Myr. Finally, the peak C objects form a broad distribution 
at faint magnitudes. 

There is one major morphological difference between the luminosity functions at these 
two bands: there is a large, persistent depression around Mk = 13-15 at all ages (Figure 
5b). This feature is located near trough F in the bottom panel of Figure 3 (the bolometric 
luminosity function); however, in Figure 5b (X-band), this feature is much deeper. Recall 
that there is a change in slope for BCk at Te// f« 1300 K (the L/T transition); this slope 
change is caused by the onset of methane formation (Figure 2). This causes the brightest T 
dwarfs to be substantially fainter at i^-band than the faintest L dwarfs. The L/T transition 
occurs over a small range of effective temperatures, leaving few objects in the transition 
region at any given age and producing the deep trough in the i^-band luminosity function. 

There is also evidence for a dip in the /-band luminosity function around Mj — 18-20 for 
all ages (Figure 5a). This feature is due to the presence of dust grains. At effective tempera- 
tures greater than ~1300 K dust exists in a homogeneous mixture in substellar atmospheres 
and smoothes out many spectral features. However, below 1300 K - the temperature also 
associated with methane formation - dust rains out of the photosphere or forms clouds (A. 
Burrows 2002, private communication; Tsuji, Ohnaka, Aoki, & Nakajima (1996)). This re- 
moves the smoothing effect, allowing the strong atomic and molecular absorption features to 
dominate. This results in a rapid decrease in flux throughout the near-infrared over a short 
range in effective temperature, which in turn creates trough F. Trough F is most pronounced 
in the X-band, because both rainout and methane formation occur around the same effective 
temperature. 

Since the effective temperature at which the L/T transition occurs varies as a function 
of mass and age, the location of the methane feature may also change as a function of 
mass and age. The bolometric correction data that we use is based on observations of field 
brown dwarfs, which are older and therefore more massive than younger objects with same 
effective temperatures. Until bolometric corrections can be obtained as a function of time 
and effective temperature it will be difficult to determine how strongly the position of the 



^We can provide luminosity functions using any desired mass function parameters for any bandpass, if 
we are given or have the appropriate BCs. Please send requests to pallenOhep . upenn . edu. 
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trough F depends on mass. 

4. Discussion 

4.1. Mciss Function Parameters 

Our models include one free parameter: a2, the power law exponent at lower masses. 
Figure 6 displays /- and M-band 100 Myr luminosity functions for three different values of 0^2 • 
As expected, the steeper mass function, corresponding to a2 = 1.2, leads to a larger relative 
number of low luminosity objects. In contrast, with a2 = 0.5, the luminosity function slopes 
significantly downward from peak A towards fainter objects. The joining mass, mi2, can 
also be varied as a free parameter. However, varying rriu throughout the range of 0.03 Mq 
to 0.1 Mq produces minimal (~5%) changes in the amplitude of the features in the output 
luminosity functions, so for the models presented here we have fixed rriu to 0.08 Mq. 

4.2. Comparison with Observational Data 

We test our models against observed cluster data for three young clusters: Upper Sco 
(Ardila, Martin, & Basri 2000), IC 2391 (Barrado y Navascues, Stauffer, & Patten 1999), and 
the Pleiades (Bouvier et al. 1998) (Table 1). Ardila, Martin, & Basri (2000) obtained RIZ 
photometry using the 60 cm Schmidt telescope at Cerro Tololo Inter- American Observatory 
(CTIO) to survey ~14 square degrees of Upper Sco, finding 138 candidates. Twenty-two 
candidates were followed up spcctroscopically and 20 were found to be cluster members. 
Barrado y Navascues, Stauffer, & Patten (1999) used the 4m Blanco and 1.5m telescopes at 
CTIO to observe ~2.5 square degrees of IC 2391 at VRIcZ to construct their 132 object sam- 
ple. They followed up their initial selections by cross referencing with the 2MASS database 
to obtain more rigorous color cuts. Bouvier et al. (1998) used the Canada- France-Hawaii 
Telescope (CFHT) to obtain 62 cluster candidates in ~2.5 square degrees of the Pleiades 
at R and /, twenty-six of which are claimed to be very low-mass stars or brown dwarfs. 
These 62 candidates were then added to an older sample of ~440 stars (Hambly, Hawkins, 
& Jameson 1993). Each group estimated their foreground/background contamination due 
to photometric selection to be 25% to 50%. Consequently, the data shown here may have 
significant contamination from non-members, especially at the faint magnitudes, and the 
conclusions we draw should be viewed as preliminary pending the availability of deeper, less 
contaminated samples. 

Figure 7 displays cluster models, with a2 — 0.5, produced for the specific ages corre- 
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sponding to each data set. Complications arise in the comparison of the data to our models 
because the brightest model magnitude bins are incomplete due to the upper mass limit of 
the Burrows et al. (2001) models. The spread in age used for each model cluster causes 
each model magnitude bin to be composed of both massive, old objects and young, less mas- 
sive objects. The complication arises at the brightest model bins because the older objects 
would have masses greater than O.ISM©, which are not included in the Burrows et al. (2001) 
models. To account for this, we discard those most luminous model bins (compare Figure 3 
with Figure 7). The brightest reliable model point is then normalized to the corresponding 
data point. The resulting match between the model predictions and the data is extremely 
encouraging, considering that a rigorous fit to the data was not performed. This gives us 
the confidence to make further predictions about future cluster observations. 

The youngest of these three clusters. Upper Sco, is at the young edge of our models. The 
fit to this cluster is quite good due to the large number of objects detected. The other two 
clusters considered, IC 2391 and the Pleiades, each have two different age estimates derived 
using two independent techniques (sec Table 1). In both cases, the younger age estimate 
(35 and 100 Myr, respectively) is based on matching the turnoff from the main-sequence 
against the appropriate stellar isochrones (such as Meynet, Mcrmilliod, & Maeder (1993)). 
The older age estimate is derived from the measurements of the lithium depletion boundary 
among low-mass cluster members. Since all of these dwarfs are fully convective primordial 
lithium is cycled through the core regions and progressively destroyed at a rate which depends 
on the central temperature. Since the latter parameter depends on mass, lithium is depleted 
more rapidly in higher mass dwarfs. Thus, as a cluster ages, the boundary between dwarfs 
with and without detectable lithium absorption moves to lower luminosities and later spectral 
types. Theoretical models can be used to calibrate that variation, and hence estimate cluster 
ages (Stauffer, Schultz, & Kirkpatrick 1998). 

We can compare the main-sequence fitting and lithium-depletion ages against our es- 
timate based on the predicted morphology of the luminosity function. In the case of the 
Pleiades, we cannot distinguish between ages of 100 Myr and 120 Myr, since the discrimi- 
nating features lie 3 to 5 magnitudes fainter than current survey limit. However, the IC 2391 
data are clearly a better match to the 35 Myr luminosity function than the 53 Myr function. 
The most notable feature in the data is the pronounced peak at Mj ~ 14 (Figure 7, middle 
panel) . This corresponds exactly to the position of peak D (high-mass brown dwarfs) in the 
35 Myr model. Peak D hes ~ 2 magnitudes fainter for the 53 Myr age favored by a recent 
lithium depletion analysis (Barrado y Navascues, Stauffer, & Patten 1999). However, the 
data displayed for IC 2391 includes both the probable and the possible members, as defined 
by Barrado y Navascues, Stauffer, & Patten (1999). The majority of the cluster stars con- 
tributing to the peak are only possible members, and further observations are required to 
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determine the extent of background contamination. 

To detect features such as peak D, the sample size for an observed cluster must exceed a 
threshold value dictated by the Poisson uncertainty in each bin. For IC 2391 this corresponds 
to ~10-15 stars per bin for a 3a difference. Therefore, since peak D covers a range of one 
and a half magnitudes (three model bins), we require there to be a total of 30-45 stars in 
the range 13.5 < Mj < 15. This is on the edge of the current data for IC 2391 with 22 stars 
in that range. 

In general, we note that the location (in magnitude) of peak D is a strong function of 
age. This may prove a sensitive, albeit model dependent, diagnostic of the age of clusters 
with ages between 15 Myr and ~200 Myr. To date, published observations of IC 2391 are 
the only data sufficiently sensitive to detect peak D. Further studies, extending to fainter 
magnitude limits in this and other clusters, are clearly warranted. 

4.3. Prospects for Future Observations 

The sensitivity limits of current young cluster studies, with the exception of IC 2391, 
preclude observations of peak D and very faint objects. The limits of the Bouvier et al. 
(1998), Barrado y Navascues et al. (2001), and Ardila, Martin, & Basri (2000) surveys are 
given by the faintest data points in Figure 7 and correspond to mj ~ 20-21. Figure 8 
shows /-band luminosity functions for clusters with a range of ages and this sensitivity. The 
hmit of detectabihty of peak D depends on the distance and age of the observed cluster. 
Additionally, peak D only exists in clusters older than 15 Myr. With current sensitivity 
limits, relatively few clusters lie close enough to probe the lower end of the luminosity 
function. Furthermore, cluster observations are often hampered by extinction (particularly at 
/-band), and determining cluster membership can be problematic. However, recent advances 
in near- and mid-infrared instruments, wide-field imagers, and space-based missions offer 
encouraging prospects to circumvent these difficulties. 

Infrared observations of very low-mass stars and brown dwarfs are advantageous because 
these objects are quite red in optical— infrared colors (I—JHKM), and incompleteness effects 
due to extinction are reduced. Additionally, at near-infrared wavelengths very low-mass stars 
and brown dwarfs are close to their blackbody peaks. Of the three primary near-infrared 
bands {JHK), J-band is best suited for these studies (Figure 9). This is because J-band 
lacks the severe methane absorption feature which causes great distortions in the /T-band 
luminosity function around the L/T transition and depresses /T-band fluxes (Nakajima et 
al. 1995). //-band is intermediate between J and K, and while the methane absorption 
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ls not as severe as that in X-band, there is some some methane absorption at iJ-band 
wavelengths. The combination of these reasons make the J-band the ideal choice for future 
cluster observations from the ground. 

Several ground-based near-infrared imagers that have recently come online provide suffi- 
cient sky coverage and sensitivity that significant portions of nearby clusters can be mapped 
efficiently. For example, Flamingos on Gemini South has a FOV = 2.6' x 2.6' and reaches 
mj = 20.3 in 300 sec (5a point source detection) (Elston 1998). WIRC-2K (Wide-field 
InfraRcd Camera) on Palomar has a FOV = 9'x9' and reaches nij = 19.9 in 300 sec (Sa 
point source detection) (Wilson et al. 2002). With these and other current and forthcoming 
instruments it will be possible to probe 3 to 5 magnitudes deeper than past surveys, detect 
peak D, and better determine the behavior of the very low-mass/substellar mass function. 

Although mid-infrared observations yield the best possibilities for detection of peak 
D and troughs B and E (Figure 10), ground-based mid- infrared observations are difficult 
because atmospheric water bands limit sensitivity. The forthcoming Space Infrared Telescope 
Facility (SIRTF) mission will probe to sufficient sensitivities to detect these intermediate 
predicted cluster luminosity function features. Figure 10 shows our M-band (4.8 fj,m) models 
with the SIRTF InfraRed Array Camera (IRAC) Channel 2 (4.5 yum) sensitivity hmits (Mm = 
13) for a cluster at 200 pc. The predicted sensitivity limit corresponds to a 5a detection of 
a niM = 19 point source in 200 seconds (Hora ct al. 2002). For this sensitivity and distance 
both troughs B and E will be easily detected by IRAC/SIRTF, as will be peaks D and C for 
clusters younger than ~200 Myr. 



5. Summciry 

We have combined theoretical evolutionary tracks for low-mass dwarfs with empirical 
bolometric correction measurements to identify several characteristic features in the pre- 
dicted luminosity functions at IJ H K M-hands. The most significant feature that we find is 
an intermediate peak in clusters that arc tens to hundreds of millions of years old (denoted as 
peak D) corresponding to high-mass brown dwarfs that have ceased to burn deuterium. This 
feature appears to be present in observations of the young cluster IC 2391, and the peak's 
location in absolute magnitude suggests an age of 35 Myr. The location of this feature may 
prove a useful chronometer for young clusters. 

Strong spectral features, such as the onset of methane absorption at the L/T transition, 
produce significant effects in the i^-band luminosity function. Both these features and peak 
D should be detectable in nearby open clusters with the current generation of wide-field 
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infrared imagers and SIRTF, allowing comparison of the model predictions with clusters 
spanning a much wider range in age. 

P.R.A. acknowledges support by a grant made under the auspices of the NASA/NSF 
NStars initiative, administered by JPL, Pasadena, CA. D.E.T. acknowledges support from 
NASA through a Space Telescope Science Institute grant to Gary Bernstein. Thanks to 

Adam Burrows for making his models available for our use and for answering our questions 
about them. P.R.A. thanks Kelle Cruz for her help in perfecting the plots in this paper. 
The authors would also like to thank the referee for his/her valuable comments and quick 
response, as well as our editor, Paula Szkody, for her timeliness and excellent choice of 
referee. 

REFERENCES 

Ardila, D., Martm, E., & Basri, G. 2000, A J, 120, 479 

Baraffe, I., Chabrier, G., AUard, F., & Hauschildt, P. H. 1998, A&A, 337, 403 

Barrado y Navascues, D., Stauffer, J. R., & Patten, B. M. 1999, ApJ, 522, L53 

Barrado y Navascues, D., Stauffer, J. R., Briceno, C. ., Patten, B., Hambly, N. C., & Adams, 
J. D. 2001, ApJS, 134, 103 

Bouvier, J., Stauffer, J. R., Martin, E. L., Barrado y Navascues, D., Wallace, B., & Bejar, 
V. J. S. 1998, A&A, 336, 490 

Burgasser, A. J. et al. 2002, ApJ, 564, 421 

Burrows, A., Hubbard, W. B., Lunine, J. I., & Liebert, J. 2001, Reviews of Modern Physics, 
73, 719 

Chabrier, G., Baraffe, I., AUard, F., & Hauschildt, P. 2000, ApJ, 542, 464 
Dahn, C. C. et al. 2002, AJ, 124, 1170 
Elston, R. 1998, Proc. SPIE, 3354, 404 
Geballe, T. R. et al. 2002, ApJ, 564, 466 

Hambly, N. C, Hawkins, M. R. S., & Jameson, R. F. 1993, A&AS, 100, 607 

Hora, J. L., et al. 2002, in Infrared Space Telescopes and Instruments, SPIE Proceedings, 
in press 



Kirkpatrick, J. D. et al. 2000, AJ, 120, 447 
Kroupa, P. 2001, MNRAS, 322, 231 
Leggett, S. K. et al. 2002, ApJ, 564, 452 

Mcynet, G., Mermilliod, J.-C, & Maeder, A. 1993, A&AS, 98, 477 
Muench, A. A., Lada, E. A., & Lada, C. J. 2000, ApJ, 533, 358 
Miller, G. E. & Scab, J. M. 1979, ApJS, 41, 513 

Nakajima, T., Oppcnhcimer, B. R., Kulkarni, S. R., Golimowski, D. A., Matthews, K., & 
Durrance, S. T. 1995, Nature, 378, 463 

Noll, K. S., Geballe, T. R., & Marley, M. S. 1997, ApJ, 489, L87 

Rcid, 1. N. & Gizis, J. E. 1997, AJ, 113, 2246 

Reid, I. N. et al. 1999, ApJ, 521, 613 

Reid, 1. N., Gizis, J. E., Kirkpatrick, J. D., & Koerner, D. W. 2001, AJ, 121, 489 

Reid, I. N. & Cruz, K. L. 2002, AJ, 123, 466 

Reid, I. N., Gizis, J. & Hawley, S., 2003, AJ124, 2721 

Salpeter, E. E. 1955, ApJ, 121, 161 

Scab, J. 1998, ASP Conf. Ser. 142: The Stellar Initial Mass Function (38th Herstmonceux 
Conference), 201 

Stauffer, J. R., Schultz, G., & Kirkpatrick, J. D. 1998, ApJ, 499, L199 

Tsuji, T., Ohnaka, K., Aoki, W., & Nakajima, T. 1996, A&A, 308, L29 

Wilson, J.C. et al., "A Wide-Field Infrared Camera for the Palomar 200- inch Telescope," in 
Proc. SPIE 4841, in press 



This preprint was prepared with the AAS IM^jX macros v5.0. 



-14- 



Table 1. Cluster Information 



Cluster 


MS turnoff age 


Li-depletion age 


Distance 


Candidates 


Ref. 


Name 


(Myr) 


(Myr) 


(pc) 


(#) 




Upper Sco 


5 




145 


138 


(1) 


IC 2391 


35 


53 


155 


132 


(2) 


Pleiades 


100 


120 


120 


62 


(3) 



References. — (1) Ardila, Martin, & Basri (2000); (2) Barrado y Navascues et al. 
(2001); (3) Bouvier et al. (1998) 
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Log Age (Gyrs) 

Fig. 1. — Burrows et al. (2001) model evolutionary tracks as a function of effective temper- 
ature and age. Solid curves show main sequence objects (M > 0.075 Mq); dashed curves 
show substellar objects massive enough to burn deuterium but not hydrogen (13 Mj < M < 
0.075 Mq); and dot-dashed curves represent objects that never experience any significant 
thermonuclear reactions (M < 13 Mj). Also marked are the ranges of MLT spectral types. 
The contours correspond to the following object masses: 1, 2, 3, 5, 6, 7, 8, 9, 10, 11, 12, 13, 
15, 16, 17, 18, 19 Mj, and 0.02, 0.025, 0.03, 0.035, 0.04, 0.05, 0.055, 0.06, 0.065, 0.07, 0.071, 
0.072, 0.073, 0.074, 0.075, 0.076, 0.077, 0.078, 0.079, 0.08, 0.085, 0.09, 0.095, 0.1, 0.11, 0.125, 
0.15 Mo (1 Mj = 1047Mo). 
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Fig. 2. — Bolometric correction curves for K- and M-bands. The filled circles show individual 
measurements. The X-band data is based on the J-band corrections in Dahn et al. (2002) 
and the M-band data is from Reid & Cruz (2002). The curves are a seventh order polynomial 
fit to the data. However, the T dwarfs (the coolest one M-band and two iiT-band points) are 
joined by straight line segments to the rest of the fit to preserve the sharpness of the L/T 
transition. The abrupt change in slope of the BC/^ data around 1300 K occurs at the L/T 
transition and is primarily due to methane absorption at temperatures cooler than 1300 K. 
For temperatures cooler than T~900 K we assume the BCs are constant. 
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Fig. 3. — Modeled bolometric magnitude cluster luminosity function for three ages 
(10, 100, and 1000 Myr) with a2 = 1.0. Each feature is labeled as discussed in §3.1. Note 
that the sharp drop on the bright side of peak A is unphysical and is a result of the high mass 
edge of the Burrows models. The normalization of the y-axis is dependent on the overall 
normalization of the distribution from which the luminosity function was created (§2.1). 
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Fig. 4. — The same as the middle panel of Figure 3 (100 Myr) for three different low-mass 
cutoffs. The black shading corresponds to objects with 0.075 Mq < M < 0.15 Mq, the 
gray shading to 13 Mj < M < 0.075 Mq, and the unshaded to M < 13 Mj. This figure 
demonstrates that the intermediate peak, D, is populated solely by objects that burned 
deuterium. 
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Fig. 5. — Cluster luminosity function evolutions for /- (a) and iC-bands (b). Light areas 
indicate relative high numbers, and dark areas indicate low numbers of objects. The same 
feature labels used in Figure 3 are used here. See the discussion in §3.2 for details on each 
feature. Each figure has the same grayscale normalization. 
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Fig. 6. — /- and M-band luminosity functions at 100 Myr for three different power laws. 
As expected, increasing 012 greatly increases the relative number of brown dwarfs compared 
to stars. For a2 = 0.5 the substellar luminosity function slopes steeply down from the main 
sequence. 
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Fig. 7. — /-band data (triangles) for Upper Scorpius (Ardila, Martin, & Basri 2000), IC 2391 
(Barrado y Navascucs ct al. 2001), and the Pleiades (Bouvier et al. 1998), compared to /-band 
models for 5, 35/53, and 100/120 Myr old clusters, respectively, with 0:2 = 0.5. Absolute 
magnitudes are calculated based on the reported distances to each cluster. We also show 
the Pleiades data from the photographic survey of Hambly, Hawkins, and Jameson 1993 
(diamonds). The models are normalized to the data points indicated by the sohd circles as 
discussed in §4.2. Also, the first two to three model bins are not plotted as discussed in §4.2. 
We show two different age estimates for IC 2391, 35 Myr (solid) and 53 Myr (dotted), and 
the Pleiades, 120 Myr (sohd) and 100 Myr (dotted). Peak D appears to be detected in the 
IC 2391 data and distinguishes between the two age estimates. 
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Fig. 8. — Model /-band cluster luminosity functions for a2 — 1.0 and ages 10, 100, and 
1000 Myr with the /-band sensitivity limit from Bouvier et al. (1998) (dotted hne). This 
limit corresponds to 90% completeness in 1200 seconds using the University of Hawaii 8K 
CCD mosiac at the CFHT for a cluster at 200 pc. Note that this sensitivity limit lies in 
the middle of the stellar/substellar trough for each of the model clusters. This implies that 
most ground-based /-band searches are unable to detect peak D except in relatively young 
nearby clusters. 
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Fig. 9. — Similar to Figure 8, but for J-band. The sensitivity limit is for a cluster at 200 pc 
using Flamingos at Gemini South which can reach mj — 20.3 in 300 sec for a 5a point 
source. 
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Fig. 10. — Similar to Figure 8, but for M-band. The sensitivity limit corresponds to a 5a 
point source detection in 200 seconds with IRAC Channel 2 on SIRTF for a cluster at 200 pc. 



